Coordinate space proton-deuteron scattering calculations including Coulomb force 

effects 
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We present a practical method to solve the proton-deuteron scattering problem at energies above 
the three-body breakup threshold, in which we treat three-body integral equations in coordinate 
space accommodating long-range proton-proton Coulomb interactions. The method is examined for 
phase shift parameters, and then applied to calculations of differential cross sections in elastic and 
breakup reactions, analyzing powers, etc. with a realistic nucleon-nucleon force and three-nucleon 
forces. Effects of the Coulomb force and the three-nucleon forces on these observables are discussed 
in comparing with experimental data. 
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I. INTRODUCTION 

Scattering observables of three-nucleon (3JV) systems, 
proton-deuteron (pd) scattering and neutron-deuteron 
(nd) scattering, are good sources of information about 
unknown aspects of the nuclear forces such as off-shell dif- 
ferences in nucleon-nucleon force (2NF) models, possible 
evidence of 3iV forces (3NFs), etc. Because of technical 
advantage of treating charged particles as beam, target, 
or detected particles over doing neutral ones, available 
data sets of the pd reaction are richer both in quantity 
and quality than those of the nd reaction. On the other 
hand, due to a mathematical difficulty in treating three- 
body systems with long-range Coulomb interactions, a 
precise calculation of the pd scattering, especially for en- 
ergies above the three-body breakup threshold (TBT), is 
one of the most challenging subjects in physics of few- 
body systems. 

In the last decade, some developments have been made 
in this problem by calculations based on the Kohn vari- 
ational principle l], 0] and on the momentum space Fad- 
deev equations [2 | using the screening and renormaliza- 



tion method 0, 



MM 



In this paper, we will present another approach to the 
pd scattering problem, which is based on integral equa- 
tions for wave functions in coordinate space. Calculations 
by this approach for non-Coulombnic 3JV systems with 
realistic 2NFs and 3NFs were performed for 3 H in Refs. 
[1, [j|, and for the nd scattering at energies above the 
TBT in Refs. [13, EH, H3] - A direct application of the 
Faddeev equation to Coulombnic 3N systems, namely 
the 3 He bound state and the pd scattering, is known to 
bring a severe singularity to the integral kernel due to the 
long-range character of the proton-proton (pp) Coulomb 
force. In Ref. [HI], Sasakawa and Sawada proposed a 
modification of the Faddeev equation to treat the sin- 
gularity by introducing auxiliary potentials that act be- 
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tween charged spectator particles and the center of mass 
(cm.) of the rest pair particles. Since we use an itera- 
tive method to solve the integral equations, in which one 
needs to operate the integral kernel on known functions 
repeatedly, it is essential to establish an accurate kernel 
operation for precise calculations. In the integral kernel 
of this Coulomb-modified Faddeev equation, which will 
be referred to as the SSF equation, the singularity due to 
the pp Coulomb potential is expected to be canceled by 
the auxiliary potentials, on conditions that three parti- 
cles are bound or no three-body breakup channel is open. 
Solutions of the SSF equation were successfully obtained 
for the 3 He bound state in Refs. [3, EH, EH and for the 
pd scattering at energies below the TBT in Refs. [17], [HI- 
However, the cancellation is not trivially expected when 
a three-body breakup channel opens. In this paper, we 
will treat this problem, and show how our approach prac- 
tically is applicable for pd scattering above the TBT. 

The next section is devoted to present notations used 
in this paper and to introduce the SSF equation for the 
pd scattering in integral equation form. In Sec. IIII1 we 
analyze a problem in the integral kernel of the SSF equa- 
tion due to the pp Coulomb force, and propose a method 
to perform numerical calculations. Then, some numeri- 
cal results for pd and nd scattering using a realistic 2NF 
and 3NFs will be presented in Sec. HVl Summary will be 
given in Sec. [V] Our iterative method [1, Ell is reviewed 
in Appendix and some useful functions and formulae 
are appeared in Appendices [B] and [Cj 



II. THREE-BODY SCATTERING EQUATION 
WITH COULOMB FORCE EFFECTS 



In this section, we will describe our notations and 
present the SSF equation by taking a proton(l)- 
proton(2)-neutron(3) system as an example. We will not 
consider spin's degrees of freedom, angular momentum 
dependence of the potentials, and 3NFs in describing 
our formalism because of simplicity. The deuteron thus 
is supposed to be a s-wave proton-neutron (pn) bound 
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FIG. 1: (Color online) The Jacobi coordinates of the three- 
body system. 



state with energy Ed. We use sets of coordinate systems 
{x k ,y k } (the Jacobi coordinates) to describe the three- 
body system defined as 



x k 



'3' 



yk = r k -^{r l + r J ), (1) 



where k) denote (1, 2, 3) or their cyclic permutations 
and Ti is the position vector of the particle i (see Fig. 
[T]). Subscripts to indicate particles will be omitted when 
there is no confusion. 

We write a three-body Hamiltonian in the cm. frame 

as 



H = H + Vi + V2 + Va 



(2) 



where Hq is the internal kinetic energy operator of the 
three-body system, 



H Q = T x (x) + T y (y) = --\7l-^ : Vl, (3) 



m 



3h2 \ 

Am v ' 



with nucleon mass m, and V k is a potential to describe 
the interaction between particles i and j consisting of a 
short-range nucleon-nucleon potential (2NP) V k s (xk) and 



the pp Coulomb potential V c (x 3 ) (= e 2 /a 
V k = V k s {x k ) + S K3 V c (x 3 ). 



(4) 



We begin with a differential form of the SSF equation 
for a three-body cm. energy E(> 0) fl3| . 



[E-Ho-Vf -u c (yi)]$i = [A*], 
[E-H -V 2 S -u c (y 2 )]$ 2 = [A*] 2 
[E-H -V 3 S -V c (x 3 )]$ 3 = [A$] 3 



(5a) 
(5b) 
(5c) 



The symbols [A$] fc in the right hand side denote 



(k = 1, 2), 



(*1 + *2) + (V C (X 3 ) - U C { Vl )) $! 

+ (y c (x 3 )- U c (2; 2 ))<i> 2 (fe = 3). 



The auxiliary potentials play different roles on each 
side of Eqs. ([5a|) - ([5c| . On the left hand side of Eqs. ([5a|) 
and (|5b[) . these potentials work to distort the spectator 
proton from a free state (see Eq. (fT4|) below). On the 
right hand side of Eq. ([5c|) . i.e. in [A$] 3 , the auxiliary 
potential is expected to cancel out the long-rangeness of 
the pp Coulomb potential V c . The latter point will be 
discussed later in this article. It should be noted that 
the auxiliary potentials in Eqs. ([5aj) - ([5c]) are eliminated 
when all the equations are summed up, which makes the 
sum $i + $2 + $3 an eigenstate of the Hamiltonian @. 

We will consider a pd scattering state of the initial pd 
momentum po, which gives the three-body cm. energy 



E = E 



E d 



where 



E n 



Zh 1 

4m 



Po- 



(8) 



(9) 



Integral form of the SSF equation, whose formal deriva- 
tion is given in Ref. [13], is 

$ fe = ~5 k ^ d {x k )F c {y k -p 0im ) + g k (E) [A$] fc , (10) 

where 6ij = 1 — Sj.j; & d (x) is the deuteron state; 
F c (y,po, rjo) is a scattering state in the Coulomb po- 
tential u c (y), which satisfies 

[T y (y)+u c (y)] F c (y;p Q , m ) = E Po F c (y;p ,r) ) (H) 

with the Coulomb parameter r/ = vipo) given by Eq. 
(IB2|) ; the operators Q k (E) are channel Green's functions 
defined as 



1 



E 



IE 



H 



V h s 



Gk(E) = I 



u c {y k ) 
(fc = l,2), 



(12) 



E 



IE 



H 



V s 

v 3 



V c (x 3 ) 
(fc = 3), 



with the parameter e being a small positive number to 
give outgoing waves. 

A partial-wave decomposition is performed by intro- 
ducing an angular function denoted as \a(x,y)), 



where $fe's are Faddeev components, and u c (y k ) is an 
auxiliary Coulomb potential acting between the particle 
k and the cm. of the pair ij, 



u c (y k ) = — {k = 1,2). 
Vk 



(6) 



\a(x,y)) = lY L (x)®Y e (y)} J ^ o 



(13) 



where L denotes the relative orbital angular momentum 
of the pair particles; £ the orbital angular momentum 
of the spectator particle; Jo(= L + I) and M the total 
angular momentum of the three particles and its third 
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component, respectively. The set of the quantum num- 
bers (L, £, J , M ) is represented by the index a. 

We define complete orthogonal sets of functions de- 
scribing the angular parts of the three-body system with 
a state index a and the radial part of the spectator par- 
ticle with momentum p and angular momentum £, 

2 F e {ri(p),py k ) (h = 1 2) 

|^fe,a) = \a(xk,Vk)) 

^PjdPVk) (k = 3), 

(14) 

where Fj(n, r) is the regular Coulomb function of Eq. 
(fBil ^HEt v(p) is the Coulomb parameter of Eq. (|B2|) . 
and jt{r) is the spherical Bessel function. The underline 
implies a dependence on the momentum p. These func- 
tions satisfy a complete relation, 



J a 



(15) 



and an orthogonal relation, 

(Fk,a\Fk,a>) = Sa,a'5{p - P ) = $a t a' , (16) 
/>oo 

where J a means / dp, and ( | ) integrations over 

~ a JO 

the variables, x and y. 

The channel Green's function is decomposed by the 
complete set to give, 



Qk(E) — / \j-k,a) Gk,L {J~k,a\ 
J a 



Here, Gul is a two-body Green's operator, 



G 



k,L 



E q + ie~ T L {x) - V k s (x) - 5 kl3 V c (x) ' 



(17) 



(18) 



where E q is the energy of the two-body subsystem given 
by 



E q = E 



and 



L ^ ^ m \ dx 2 x dx 



3h 7 h 2 
4m m 



2 d L(L+1) 



(19) 



(20) 



III. COULOMB FORCE EFFECTS IN THE SSF 
INTEGRAL KERNEL 

The SSF integral equation presented in the previous 
section has a form of the inhomogeneous linear equa- 
tion. We are going to solve this by applying an itera- 
tive method developed in Refs. d, [l^ ( an d references 
therein), which is called as the Method of Continued 
Fractions (MCF). In general, iterative methods to solve 



an linear equation require to operate the kernel to func- 
tions that are given in preceding iterative steps. The 
MCF algorithm, which is reviewed in Appendix lAl also 
includes such operations as indicated in Eqs. (|A8|) and 
(|A9[) . Calculations of the SSF integral kernel consist of 
two parts: a particle exchange operation and the opera- 
tion of the Green's functions. Some technical notes of the 
former part are given in Refs. [22j, |23j], and those of the 
latter part for the nd scattering above the TBT in Ref. 
[To| . which are useful also in the pd scattering. In this 
section, we will study some problems of the SSF integral 
kernel proper to the pd scattering problem. 



A. SSF integral kernel 

Let us consider to operate the SSF kernel on given 
functions Q k (x,y): 



Q k {x kl y k ) = g k {E) [A$] A 



(21) 



The channel Green's function for k = 1 or 2, where 
the pair is a pn system, possesses a pole corresponding 
to the deuteron bound state. In order to treat this pole, 
we apply a standard subtraction method, in which we 
use an identity, 



l = £> o d )(0%)| 



1-^|« O ^)(0%)| 



, (22) 



where (f> d (x) is the radial part of the deuteron wave func- 
tion with orbital angular momentum Lq(= 0), and the 
index ctQ — (Lq, £o, Jo, Mq) denotes the three-body par- 
tial wave states that couple to the two-body state with 
Lq. By applying the identity to Q k , we obtain 

G k (E) = ^|a ^)GcVo(£ P o)(<Ao| 



1 



Eq — Ed 



(k = l,2). 



(23) 



Here, Gc,t (E Po ) is the partial wave component of the 
Coulomb Green's function for the outgoing proton, 



Gc,i (E p 



1 



E P0 +ie-T eo (y)~u c (y) 



(24) 



with 



37T 



2 d £{£+!) 



4m \dy 2 y dy //- 



(25) 



The function <d k (x,y) thereby can be written as 



4 



' £M0 d (^ o ( y ) 

J a_ 

(fc=l,2), 

|-^fe,a) O k ,a(x) (k = 3). 



(26) 



Here, the function fi k ^ ag {y) (k — 1,2) represents an 
elastic component in the scattering, 



Jo 



with 



G c ^(y,2/';i? P ) = (y|Gc,K^)l2/'), 



(28) 



and a source function wjj. 6 ^ (y) given by 

41 (y) = <<Ao|[A$] fe ) 

= (0%>|T4 S | $ i + * J ->. 



(29) 



The explicit expression of the Green's function 
Gc,e(y,y';E p ), Eq. (|B5|) [20[ | . gives the asymptotic form 



(e) ^ e^o^)^fa, P0 ,) (e) 



(30) 



where Tfc2 is an amplitude defined by 



ji(e) _ _ 



4m \ f 00 ( e ) 
-T2— / dyF e {r]o,poy)yuj { k e> (y). (31) 
3ft po/ Jo 



Above the TBT, the source function 0J ka (y) reveals a 

long-range behavior of 0{y~ 5 / 2 ) even in the case of the 
nd scattering due to the particle exchange with breakup 
channel. This property was studied to develop a numer- 
ical treatment in Ref . [l(| • 

The coefficient Cfc )So (k — 1,2) and the function 
0k,a{x) (k = 1, 2, 3) in Eq. (|26|) are defined as follows: 



E q — Ed 



(<t> d rk& \v£\*i + * j ) 



h,a(x) = (x\G k ,L\Uk,a), 



(32) 



(33) 



where a source function ujk, a {x) is composed of a contri- 
bution from the short-range potential and one from the 
Coulomb potentials, 



with 



Uk,a(xk) = UJ ka (x k ) + S k ,3LOai x k) 



s h Jx k ) = (T k ^\V k s \^ + <P 3 ) 



(34) 



= V k s (xk)(FkJ$i + $ j ), (35) 

"*(X3) = (^3, a |^ C (z3)-« C (yi)|$l) 

+ (^|V c (a;a)-u (ift)|$2>. (36) 
Note that the apparent singularity of Ck.a for E q 
/.',;. orp= (E+\E d \), (see Eq. (JTSJl) is canceled 

by that of the function 8k, a {x) arising from the two- 
body Green's function, G k> L, and therefore the standard 
quadrature can be applied to perform the p-integration 
of Eq. (|2"6")) as far as both terms are treated together as 
demonstrated in Ref. [l(|. 

In actual calculations of the functions 0k,a(x), we con- 
sider an ordinary differential equation that is transformed 
from Eq. (33J, 

[E q - T L {x) ~ V k s {x) - 5 k ^V c {x)] 9 k ,a(x) = u k , z (x). 

(37) 

A boundary condition to get physical solution of this 
equation depends on energy of the two-body sub-system 
E q , and thus on the integral variable p in Eq. f|26[) via 
Eq. (fl9|) . According to the sign of E q , the range of p 
(0 < p < oo) is divided into two regions: (i) < p < 

p c = J |pr-E, where E q > 0, and (ii) p c < p < oo, where 



E q < 0. Corresponding boundary conditions for k 
are: 



0k,a(x) 



h { L + \qx) (0<p< Pc ), 
ht(i\q\x) {pc<P < oo), 



1,2 



(38) 



where h^(r) is the spherical Hankel function with the 
outgoing wave. For k = 3, where the pp Coulomb poten- 
tial is acting, we have 



«i (7(9) .9^) 



(0<P<Pc), 



6 3 ,a(x) OC { (39) 

^M^Bm ( Pc <p< 00), 

where 7(g) is given by Eq. (|C6|) and W KtP (z) is the Whit- 
taker function [2lj]. We solve Eq. (f37|) with above con- 
ditions by applying usual techniques as in the two-body 
problem, e.g. the Numerov algorithm [23| . Treatments 
of Eq. ([37)) for k = 1, 2 in the region (i), which are same 
as for the nd scattering, are described in the Appendix 
B of Ref. jlCj |. While those for the k = 3 case, where 
we need to consider Coulomb force effects, arc given in 
Appendix [U] of this article. 

The asymptotic form of the function Q k (x,y) is ob- 
tained by evaluating Eq. (|2"6"|) with the saddle-point ap- 
proximation 0, [25[ together with an explicit asymptotic 
form of 03,q(x) for < p < p c given by Eq. (|C26|I . We 
notice that the Coulomb force effects appear in the spec- 
tator variable y k for k = 1, 2 and in the pair coordinate 
ajfc for k = 3. The result is 



-e*' 



\a)i 



4K Q 



3/2 



5 



o j(iC H-i5fc,37(«) ln(2gx)-5fe, 3 »7(p) ln(2p»)) 



xB fc ,a(e), 



i? 5 /2 



(40) 



where the limit is considered to be taken with x/y being 
fixed, a hyper radius R and a hyper angle O are intro- 
duced as 



R 



x 2 + ^y 2 , 



x = i?cosO, y 



i?sinO, 



Kq and the momenta, q and p, are given by 



q = K cos 6, 



P : 



-Kq sin 0. 



(41) 



(42) 



(43) 



(44) 



Here, -Bfc, Q (0) is a breakup amplitude defined as 
1 m 1 



B k , a (Q) = - 



p h 1 - iK, L {q) 



(^kA<l)\Uk,a), (45) 



where ipk,L{x;q) is a two-body scattering solution with 
the standing wave boundary condition and K,l{q) is a 
scattering if -matrix for the two-body scattering (see Ap- 
pendix [Cj . 



B. Coulomb long range effects 

In solving the differential equation (|37[) numerically, 
we need to set a value xm by a condition that the 
source function ujk,a{x) should vanish so that the solution 
reaches its asymptotic form given by Eqs. (f3"5)) or ((551) 
for x > xm- The range of tOk,a(x) thus is an important 
issue in our calculations. Eq. (|35| shows that the range 
of the short-range potential term wj? a {x) for k = 1,2,3 

is determined by the range of V^{x). Therefore, in the 
case of k — 1, 2, where there is no contribution from the 
Coulomb term, we set xm to be a value larger than the 
range of the 2NP, e.g. 10 fm. 

In the case of k = 3 on the other hand, the source 
function includes the Coulomb term 0^(2:3), whose range 
depends on a factor 



{y c (x 3 )-u c (yi)}<i>i( : r 1 ,y 1 ) 
1 1 

%3 2/1 , 



+ (1^2). 



+ (1^2). 



In our iterative scheme (see Appendix |A"|) . the zeroth 
order of the source function uj^ a (x 3 ) is calculated by 
putting $i(a?i,yi) = <& d (xi)F c (y^ p , 770), in which the 




V C ( X3 ) X 3 



FIG. 2: (Color online) Jacobi coordinates (xi, yi, and x$) to 
describe the pp Coulomb potential V c (x^) and the auxiliary 
Coulomb potential u c \yi). 



magnitude of the variable Xi is restricted within the 
range of the deuteron size. Using an expression given 
by the definition of the Jacobi coordinates, Eq. ([1]), (see 
also Fig. 02), 



we can easily show that 

111 1 

X3 



-xi, 



1 

yi 



1 

X3 



\x 3 



Xi\ x 3 - 



(46) 



0{xf). (47) 



The same situation holds for the replacement of (1 <-> 2). 
The source term u)^(x 3 ) therefore supposed to be a short- 
range function due to a cancellation between V c and u . 

An example of the cancellation is shown in Fig. [3] (a) , 
where we plot components of ^o^ a {x 3 ) for a partial wave 
state of 1 So(pp) — Sx/2 and the total three-body angular 
momentum and parity of l/2 + , and p = 0.30 fm -1 for 
the pd scattering at incident proton energy E p = 13.0 
MeV using the Argonne Vis (AV18) 2NP [2f3]. In the 
figure, a component due to the 2NP, a (xs) (the solid 
curve), and components due to the Coulomb potentials, 
w a (X3), including only V c (the dotted curve) and both 
of V c and u c (the dashed curve) are plotted. For £3 < 2 
fm, only the term W3 a {x-s) is plotted since the Coulomb 
contributions are very small in this region. As shown by 
the dotted and dashed curves in the figure, the contribu- 
tion of V c is well canceled by that of u c for large values 
of X3 in the zeroth order calculation. 

On the other hand, Fig. [3] (b) shows the components 

of uil^ a (x3) calculated from functions &^\xk, %)k) that 
are obtained by the operation of the kernel to the ini- 
tial state. Once the integral kernel is operated, the re- 
sulting functions include the three-body breakup compo- 
nent as expressed by Eq. (f4*0"|). and thus, the range of 
x\ in such functions is not restricted to the range of the 
deuteron. As a result, the cancellation as Eq. (|4"T)l is 
no more expected for higher order calculations. This is 
demonstrated by the fact that dashed curve in the figure, 
which denotes the source term w^(a; 3 ) including both of 
V c and u c contributions, remains to be non-negligible 
for a large value of x$ . 
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(a) 



t 









K (fm) 



x 3 (fm) 



FIG. 3: (Color online) Examples of the source functions (a) 
for the zeroth order and (b) for the first order calculations. 
The solid curves denote the contribution from the short-range 
2NP u>i l(x (x3), the dotted curves that from the Coulomb po- 
tentials u>a(x3) neglecting the contribution from u c , and the 
dashed curves (2:3) including the contribution from both of 
V c and vF . Note that the scales of the vertical axis change 
at X3 = 2.0 fm. 



To include the long-range effect of ^(13) as much as 
possible, one needs to increase the value of xm much 
larger than 10 fm, which is the standard value in the 
nd calculation. This makes pd calculations much harder 
than the nd calculations. Since Fig. [3] implies that the 
source function u>3 j0l (x3) is dominated by the 2NP con- 
tribution uj§ a (x3), we decide to include the effect of the 
long-range contribution partially by multiplying oj^^x^) 
by a cutoff factor 



(xs/RcV 



(48) 



for higher order than the zeroth order in our iterative 
procedure of the MCF. 

Validity of this procedure is examined in TablelH where 
we compare some results of the pd eigenphase shift and 
mixing angle parameters in convention defined in Ref. 
[27l ] . for a partial wave state with the total angular mo- 
mentum and parity of l/2 + calculated with the AV18 
potential. In the table, the last column shows the results 
by the Kohn variational principle (KVP) 28|. The rest 
columns do our calculations. The first column denoted 
as WCO shows results calculated by completely neglect- 
ing ujC(x 3 ) in all order calculations of the MCF iteration 
scheme to solve the SSF equation. Calculations with 
keeping oj^(xs) only for the zeroth order without cut- 
off but neglecting ui^(x 3 ) for higher order calculations 
are shown in the second column (WCn). In the third to 
fifth columns, we show results with the cutoff for non- 
zeroth order by choosing N = 4 and Rc = 4 fm, 6 fm, 
and 8 fm, respectively. A comparison of the numbers in 
the table indicates that effects due to the neglect of the 
long-range term uj^ (£3) in the integral kernel may be an 
order of a few % in the phase shift parameters, and sug- 
gests that the partial inclusion with the cutoff factor with 
(N,Rc) = (4, 8 fm), e.g., may produce sufficient results. 
We remark that the imaginary part of the 4 S]y2 parame- 
ter reveals a rather slow convergence, which might affect 
scattering observables. This point will be discussed later. 



IV. NUMERICAL RESULTS 

In this section, we will present some numerical results 
by the formulation described in the previous sections. 
Technical details of introducing spins' degrees of free- 
dom, 3NFs, etc., are given in Refs. @, E3, Hi SI • As a 
standard 2NF model, we choose the AV18 potential [26| . 
Three-nucleon partial wave states, which the 2NF and 
3NFs act, are restricted to those with total two-nucleon 
angular momenta J < 6 for bound state calculations 
and J < 4 for scattering calculations. In scattering cal- 
culations, total 3N angular momentum is truncated at 
Jo = 19/2, while 3NF's are switched off for 3N states 
with J > 13/2. 

As described in the previous section, the Coulomb 
source term w^(x) in the SSF integral kernel is treated 
by multiplying with the cutoff factor Eq. (|48[) for higher 
order in the MCF iteration. Comparisons of calculations 
performed by taking three different sets of (N, Rc) in 
Table U show that a satisfactory convergence is obtained 
with parameters of (N, Rc) = (4, 8 fm) for elastic ob- 
servables, and we thus proceed with these parameters 
referring to them simply as pd calculations. A conver- 
gence problem for three-body breakup observables will 
be discussed in a subsection below. 

Calculated binding energy of 3 H ( 3 He) with the AV18 
potential is 7.626 MeV (6.928 MeV), which is under- 
bound by about 1 MeV compared to the empirical value 
of 8.482 MeV (7.718 MeV). It is well known that a 3NF 
that caused by the exchange of two pions among three 
nucleons (27rE-3NF) produces enough attraction to ex- 
plain the empirical binding energy. In this paper, we use 
a new version of the Brazil 27rE-3NF [ll| with a dipole 

' A 2 -rr 



form factor of the cutoff mass parameter A, ■ A 2 +q 2 

for the ttNN vertex (BRa). In a combination with the 
AV18 2NP, A is chosen to be 660 MeV (AV18+BR 660 ) to 
give the binding energy 8.492 MeV (7.763 MeV) for 3 H 
( 3 He). 



A. Differential cross section in elastic scattering 

First, we compare calculations approximately includ- 
ing the pp Coulomb force effects with those of the pd 
calculations for the differential cross section a(9) of the 
pd elastic scattering, where 9 is the scattering angle in 
the cm. system. We take two approximate calculations: 
one is the WCO calculation, which is presented in the 
previous section. The other one, which will be denoted 
as APn, is an approximate calculation, in which the scat- 
tering amplitude due to the short-range 2NF is replaced 
by a corresponding nd scattering amplitude [29(. It is 
expected that the WCO calculations are better approxi- 
mation at lower energies since breakup effects are smaller. 
On the other hand, the APn calculations are expected to 
be better for higher energies. Fig. 2J where WCO (bold 
curves) and APn calculations (thin curves) of differential 
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TABLE I: The real and complex part of the pd eigenphase shift and mixing angle parameters, which are given in degrees, for 
J n — l/2 + state with the AV18 potential. See the text for the meaning of the calculations. 
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4 D 1/2 


(-5.33, 0.01) 


(-5.44, -0.01) 


(-5.44, -0.01) 


(-5.45, 0.00) 


(-5.45, 0.00) 


(-5.43, 0.004) 


4 Sl/2 


(-42.4, 1.93) 


(-42.4, 1.93) 


(-41.8, 2.36) 


(-41.8, 2.30) 


(-41.8, 2.15) 


(-41.8, 1.74) 


Vl/2+ 


(0.97, -0.04) 


(1.01, -0.03) 


(1.05, -0.03) 


(1.05, -0.04) 


(1.05, -0.04) 


(1.05, -0.03) 


E p = 10.0 MeV 














4 D 1/2 


(-7.15, 0.24) 


(-7.32, 0.21) 


(-7.32, 0.22) 


(-7.33, 0.22) 


(-7.34, 0.22) 


(-7.30, 0.24) 


4 Sl/2 


(-61.3, 11.6) 


(-61.5, 11.5) 


(-61.0, 12.4) 


(-60.9, 11.9) 


(-60.8, 11.9) 


(-60.6, 11.7) 


Vl/2+ 


(0.96, 0.03) 


(0.98, 0.05) 


(1.02, 0.05) 


(1.02, 0.04) 


(1.01, 0.04) 


(1.01, 0.06) 



Ref. [21 




between the data [30j, |3l|, [33] and calculations denned by 
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e (deg) 

FIG. 4: (Color online) Differential cross section of pd elastic 
scattering normalized by the pd calculation with the AV18 
potential. The bold curves represent the AV18-WC0 calcula- 
tions, and the thin curves the AV18-APn calculations. The 
solid curves denote calculations at E p = 5.0 MeV, the dashed 
curves E v = 10.0 MeV, and the dotted curves E p = 28.0 MeV. 



cross sections normalized by those with the pd calcula- 
tions at E p = 5.0, 10.0, and 28.0 MeV are plotted, looks 
to exhibit roughly these tendencies. It is remarkable that 
deviations of the APn calculations are rather large, about 
10 % even at backward angles. 

In Fig. pd calculations for a(9) of the pd elastic 
scattering at E p = 5.0 MeV, 10.0 MeV, and 28.0 MeV 
are plotted together with experimental data [3(1 l3l| . 
Comparing the calculations with the AV18 potential (the 
dashed curves) to the experimental data at 8 ~ 120°, 
where <j{9) takes the minimum, one finds that the calcu- 
lations overestimate the data at lower energies and un- 
derestimate at higher energies. The introduction of the 
27rE-3NF as shown by the solid curves reduces almost all 
of the discrepancies at lower energies. This systematic 
difference between the 2NF calculations and the data, 
which is referred to as "Sagara discrepancy" , was pointed 
out in Ref. [3(| using the APn calculations. To study 
this discrepancy in detail, we plot a relative discrepancy 







a expl 



(49) 



in Fig. [5J where 8 m in is the scattering angle where the dif- 
ferential cross section takes the minimum. The pd calcu- 
lation shows that a systematic discrepancy still remains 
when effects of the Coulomb force are treated properly, 
but with shifting transition energy from the overestima- 
tion to the underestimation to a higher energy of about 
E p = 20 MeV as compared to that by the APn calcula- 
tion, about 5 MeV. This tendency is consistent with the 
results reported in Refs. 0, 0| . 



B. Phenomenological three- nucleon force 

In Ref. [18j, it is pointed out that the introduction 
of the 27rE-3NF causes an undesirable effect to the ten- 
sor analyzing power T2i($) of the pd elastic scattering 
at energies below the TBT. Also, the 27rE-3NF is known 
to give little effect on the vector analyzing power A y (8), 
for which there exists rather large discrepancy between 
experimental data and calculations ( ll A y puzzle"). These 
facts, which are also demonstrated in Fig. [71 suggest that 
the 27rE-3NF is insufficient to comprise a Nuclear Hamil- 
tonian in addition to the realistic 2NF. Since no possible 
mechanism to produce additional 3NF to remedy above 
defects is established, at the moment, a phenomenologi- 
cal 3NF model is introduced [13, which has a form that 
typical components in 2NP: central, tensor, and spin- 
orbit components, are modified in the presence of third 
nucleon. The explicit form of the 3NF is 



yphc 



r G . 



V + V T S T (ij)Pi 



^ s e-^£[V(Si + Si)]A 



(50) 



i<j 



where Sriij) is the tensor operator acting between the 
nucleon pair ij, P\\ is the projection operator to the 



8 




60 120 

e (deg) 



180 



FIG. 5: (Color online) Differential cross sections of pd elastic 
scattering at E p = 5.0, 10.0, and 28.0 MeV. The dashed curves 
denote calculations with the AV18 potential and the solid 
curves those with the AV18+BR,660 potential. Experimental 



data are from Ref. 

Ref. or 



301 for E v 



5.0 MeV and 10.0 MeV, and 



for E p = 28.0 MeV. 
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FIG. 7: (Color online) (a) Proton vector analyzing power 
A y (8) and (b) deuteron tensor analyzing power T2i(9) of pd 
elastic scattering at E p — 3.0 MeV. The dotted curves denote 
calculations with the AV18 potential, the dashed curves those 
with the AV18+BR660 potential, and the solid curves those 
with the AV18+BR 8 oo+V p '" i potential. Experimental data 
are from Refs. [30], H- 



by about 1 MeV. The parameters are decided to repro- 
duce the following observables: the triton binding energy, 
vector analyzing power A y (9) and tensor analyzing power 
T 2 i(0) in pd scattering at E p = 3.0 MeV. The results are 
V = 25 MeV, V T = -40 MeV, and V ls = -16 MeV. 

In the present work, we will use a new version of Brazil 
3NF [ll| , which is more attractive in the 3N bound states 
than the earlier version (33[. We thus need to retune 
the value of Vo to be 36 MeV, but without changing the 
values of Vr and Vu- Calculated binding energy with 
this set of potentials (AV18+BR 800 +V ? '' ie ) is 8.482 MeV 
(7.757 MeV) for 3 H ( 3 He), and the results for A y {9) and 
T2i(9) at E p = 3.0 MeV is shown by solid curves in Fig. 



m 



C. Polarization observables in elastic scattering 



FIG. 6: (Color online) Discrepancy of the pd differential cross 
section minimum A m in defined in Eq. (|49|l using the exper- 
imental data [301 l3ll I32I ]. The solid squares denote the pd 
calculations with the AV18 potential, the triangles those with 
the AV18 + BRe60 potential, and the crosses the APn calcu- 
lations with the AV18 potential. 



spin and isospin triplet state of the pair ij, and p — 
I ( r i2 + r 23 + r li)- The range parameter tq was taken 
to be 1.0 fm, and a/ s to be 1.5 fm _1 . In Ref. [T3], the 
strength parameters, Vo, Vr, and Vi s , are determined in 
the following manner: We choose the combination of the 
AV18 potential for 2NF and a former version of Brazil 
model H| with A = 800 MeV for the 2ttE-3NF as a 
starting interaction, which makes the triton overbound 



In Ref. [I2J , it is shown that the use of the phenomeno- 
logical 3NF together with the AV18+BR 800 , which is 
tuned to reproduce the 3N binding energy; A y (9) and 
T2i(9) at E p — 3.0 MeV, is also successful in describing 
the neutron vector analyzing power A y (9) of the nd scat- 
tering at higher energies. In Fig. [51 calculations of A y (9) 
of the pd and nd scattering with the AV18+BR 80 o+V ?, ' le 
potentials are compared with experimental data at some 
energies above the TBT [13, H [Hi HI] ■ While the calcu- 
lations of the nd-A y (9) agree with the experimental data 
in similar manner as in Ref. those of the pd-A y (9) 
overestimate the data at the maximum region 9 ~ 130° 
as the energy increases. In another aspect, the calculated 
difference between the nd- and the pd-A y {9) at the max- 
imum region is decreasing as the energy increases, which 
is contradictory to the tendency of the experimental data. 
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FIG. 8: (Color online) Nucleon vector analyzing power A y (ff) of pd and nd elastic scattering at En = (a) 5.0 MeV, (b) 10.0 MeV, 
(c) 14.0 MeV, and (d) 16.0 MeV. Solid (dashed) curves denote pd (nd) calculations with the AV18+BR 8 oo+V p ' le potential. 
Experimental data are from Refs. [30l. |34| for pd (solid squares) and Refs. [351. l3q| for nd (open circles). 



In Fig. [9l calculations of the deuteron tensor ana- 
lyzing power T 2 i(#) of the pd scattering at E p = 10.0 
MeV and 28.0 MeV are compared with experimental data 
[3ll H3] ■ As in the case at the low energy, the introduc- 
tion of the 27rE-3NF shifts the calculations to wrong di- 
rection from the experimental data around 9 = 90°, and 
the phenomenological 3NF works to reproduce the data. 
Another interesting feature appears at 9 ~ 130°, where 
T2i(0) takes the maximum as follows: while the calcu- 
lations with the AV18 potential deviate from the data, 
the introduction of the both 3NFs results the almost the 
same and remedy the discrepancy well. These suggest 
that Ti\(9 ~ 90°) is sensitive to tensor components of 
nuclear forces and T 2 \{9 ~ 130°) to central components. 

In Fig. [T0l energy dependence of the deuteron ten- 
sor analyzing power T^iiff) at 9 = 90° in the pd elastic 
scattering for calculations with the AV18, AV18+BR 6 6o, 
and AV18+BR 800 +V p ' l ^potentials is shown m compar- 
ing with available data [30, HH [U, [3?J • The figure shows 
that the introduction of the phenomenological 3NF is still 
consistent with data at higher energies. It is interesting 
if further T 21 (9) data at E p = 20 MeV to 30 MeV, where 
experimental data are missing, are consistent with the 
calculation or not. 

Polarization-transfer coefficients are another interest- 
ing observables, which are sensitive to spin-dependent 
interactions. In Fig. QTJ the polarization-transfer coef- 
ficient Ry (9) of pd and nd elastic scattering at Em = 
19.0 MeV are compared with experimental data [38l l39l| . 
One interesting point, which has been already remarked 



in Ref. [4(j, is that the Coulomb force effect in the cal- 
culation is opposite to that in the data at 9 ~ 110°. In 
addition, the figure shows that the experimental data in- 
dicates that the AVI8+BR660 potential is favored than 
the AV18+BR 80 o+V p/le , implying this observable may 
be useful in distinguishing various 3NF models that re- 
produce other observables equally. 

D. Breakup cross section 

Finally, we will show some results for differential cross 
sections of kinematically complete three-body breakup 
reactions, d(p,pp)n and d(n,nn)p^ which are character- 
ized by configurations of three particles in the final state. 
Here, we will discuss four different kinematical conditions 
that include the following typical configurations, whose 
experimental data at En = 13.0 MeV are available for 
the pd-breakup in Ref. .41] and for the nd-breakup in 
Refs. dHH: 

(a) Collinear (COL) configuration, in which three nu- 
cleons align on a line with the unobserved nucleon being 
at rest in the cm. system; 

(b) Final state interaction (FSI) configuration, in 
which the relative energy between the unobserved nu- 
cleon and one of the observed nucleon is zero; 

(c) Space star (SST) configuration, in which three nu- 
cleons have equal energies and interparticle angles of 120° 
in the cm. system, and the plane spanned by the three 
nucleons is orthogonal to the beam axis; 



10 



0.05 


(a) 










eg 
| — 




k \ 


0.00 






-0.05 




it 

V 

V-'/f 



e (deg) 



0.2 



-0.2 



.(b) 








V P 

\ L : 





60 , , 120 180 

e (deg) 

FIG. 9: (Color online) Deuteron tensor analyzing power 
721(6*) of pd elastic scattering at E p = (a) 10.0 MeV and (b) 
28.0 MeV (or equivalently E d = 20.0 MeV and 56.0 MeV, 
respectively, in a deuteron incident scattering). The dot- 
ted curves denote calculations with the AV18 potential, the 
dashed curves those with the AV18+BR660 potential, and the 
solid curves those with the AvT8+BR 8 oo+V p ' le potential. Ex- 
perimental data are from Ref. [37| f° r E p ~ 10-0 MeV and 
Ref. [H for E p = 28.0 MeV. 



(d) Quasi-free scattering (QFS) configuration, in which 
the unobserved nucleon is at rest in the laboratory sys- 
tem. 

First, we have checked a convergence of the breakup 
cross sections with respect to the cutoff procedure of the 
long-range Coulomb force effect with Eq. (|48|) . Calcula- 
tions with three parameter sets shown in Table [H namely 
(N,R C ) = (4, 4 fm), (4, 6 fm), and (4, 8 fm), for the 
SST and the QFS configurations agree with one another 
excellently, however those for the COL and the FSI con- 
figurations do in part as shown in Figs. [121(a) and (b). 

The visible deviations in Figs. [12] (a) and (b) appear 
at a kinematical condition where the relative pn energy 
is small, i.e., E pn < 0.5 MeV, which might be caused by 
a small change in the pp interaction due to our Coulomb 
treatment. To check this peculiar behavior, we have in- 
vestigated the dependence of the pn-FSI cross sections 
in the nd breakup reaction on the neutron-neutron (nn) 
interaction using two different 2NP models: a charge in- 
dependent 2NP, Argonne V i4 (AV14) 0, in which the 
nn force is equal to the pn force in the 1 Sq state, and 




10 20 30 

E p (MeV) 

FIG. 10: (Color online) Energy dependence of deuteron ten- 
sor analyzing power T2i(8 — 90°) of pd elastic scattering. 
The dotted curve denotes calculation with the AV18 poten- 
tial, the dashed curve one with the AV18+BR660 potential, 
and the solid curve one with the AV18+BR 800 +V p ' le poten- 
tial. Experimental data are from Ref. [3J] (squares), Ref. 
[13] (circles), and Ref. [HI (triangle). 
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FIG. 11: (Color online) Polarization-transfer coefficient 
Kjj (9) of pd and nd elastic scattering at Em — 19.0 
MeV. The dotted curve denotes the pd calculation with 
the AV18 potential, the dashed curve the pd one with the 
AV18+BR660 potential, the solid curve the pd one with the 
AV18+BR 8 oo+V p,l ' ! potential, and the bold curve the nd one 
with the AV18+BR 80 o+V p ' le potential. Experimental data 
are from Ref. [38| ] for pd (solid squares) and from Ref. j39( ] 
for nd scattering (open circles). 

its modified version (AV14') made in Refs. [H EH by 
considering an charge-dependent potential to distinguish 
the nn force from the pn force. Results are shown in 
Figs. [12] (c) and (d), which demonstrate that a change 
in a nn force actually results in non-negligible effects for 
the pn-FSl cross sections. 

Next, as a reference, cross sections of the FSI configu- 
ration for d(p,pn)p and d(n, np)n reactions at En = 13.0 
MeV are plotted in Fig. [T^l In this configuration, the pn- 
FSl occurs around S — 3 MeV and the pp-FSI or rm-FSI 
does around 5=11 MeV. The pp-FSI cross sections are 
suppressed by the pp Coulomb force compared to the nn- 
FSI cross sections, but not completely. This may suggest 
that we need to improve the Coulomb cutoff procedure 
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possibly by extending the range of the cutoff function to 
treat the pp-FSI more correctly. 

From these considerations, we conclude that our cal- 
culations successfully converge for the most of breakup 
configurations possibly except for a limited region with 
relative energy of two nucleons being close to zero. 

In Fig. 1141 results of the pd- and nd-breakup cross 
sections for the above four configurations with the AV18 
and the AV18+BR,66o potentials are compared with the 
experimental data. Effects of the Coulomb force are vis- 
ible for the COL, SST, and QFS configurations, but not 
so for the FSI configuration, which is consistent with the 
result of the momentum space calculations @ . On the 
other hand, effects of the 27rE-3NF are small except for 
the QFS configurations. In the momentum space ap- 
proach ((J, three-nucleon force effects are incorporated 
alternatively in terms of an explicit introduction of a sin- 
gle virtual A-isobar excitation. Their results also show 
that effects of the A-isobar in the breakup cross sections 
are small for the COL, FSI, and SST configurations, and 
are visible for the QFS configuration. 



APPENDIX A: METHOD OF CONTINUED 
FRACTION 

In this appendix, we summarize the MCF algorithm 
to solve the SSF equation (see Refs. [1, [l^] and refer- 
ences therein). Let us consider to solve a linear integral 
equation 



|$) = \F)+QA\$) 



(Al) 



In the notation of the present work, |<E>) and \F) are 
expressed as vectors, 



' $(a:i,tfi) 
*(^2,y 2 ) 

\ $(^3,2/3) 



^ d (x 1 )F c (y 1 ;p ,7 l0 ) 
\F) = | $ d (x 2 )F c (y 2 ;p ,r, ) 




(A2) 



(A3) 



V. SUMMARY 

We presented a practical method to solve the pd scat- 
tering problem at energies above the threshold of the 
deuteron breakup in accommodating effects of the long- 
range pp Coulomb force as accurate as possible. Al- 
though the convergence with respect to the cutoff pro- 
cedure of the long-range Coulomb force effect is left as a 
future problem at a particular kincmatical condition of 
breakup reactions, a successful convergence is obtained 
for elastic observables and for the most of kincmatical re- 
gions in breakup reactions. We thereby calculated some 
observables in pd and nd reactions at energies up to 30 
MeV. Effects of the two-pion exchange 3NF and the phe- 
nomenological 3NF to reproduce low-energy 3N observ- 
ables are examined for pd observables at higher energies, 
and then some discrepancies between calculations and 
experimental data as well as inconsistencies between cal- 
culations and data with respect to Coulomb force effects 
are observed. Studies for searching realistic mechanisms 
to produce interaction models to remedy these defects 
with calculations by the formalism presented in this pa- 
per are now in progress. 



and G and A are as matrices 
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Setting and as 

F^) = |F), 

g m = g, 

we define \F®) and (i = 1, 2, . . . ) as follows 

: gw 



F 1 ' 



(A6) 
(A7) 

(A8) 



A 1 ( 



F 



[o] 



Acknowledgments 

The author thanks H. Paetz gen Schieck for providing 
experimental data. Numerical calculations in this article 
are supported by Research Center for Computing and 
Multimedia Studies, Hosei University. 



EI^)7FHRWliT( F 1- (A9 » 



A\Fti-y) 

Introducing |<i>M) (i = 0, 1, 2, . . . ) as solutions of 



$W = 



g^A 



(AlO) 



12 




S (MeV) S (MeV) 

FIG. 12: (Color online) Differential cross sections of pd breakup reactions (a) for the COL configuration and (b) for the FSI 
configuration; those of nd breakup reactions (c) for the COL configuration and (d) for the FSI configuration, at En = 13.0 MeV. 
In (a) and (b), dotted curves denote calculations with (N,Rc) = (4, 4 fm) for the AV18 potential, dashed curves (N,Rc) = 
(4, 6 fm), and solid curves (N,Rc) = (4, 8 fm). In (c) and (d), solid curves denote the calculations with the AV14 potential 
and dashed curves the AV14' potential. 
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FIG. 13: (Color online) Differential cross sections of d(p,pn)p 
and d(n, np)n reactions for the FSI configuration at En = 
13.0 MeV with the AV18 potential. Dotted curves denote 
calculations for the pd breakup reaction with (N,Rc) = (4, 
4 fm), dashed curves (N,Rc) = (4, 6 fm), and solid curves 
(N,Rc) = (4, 8 fm). Dot-dashed curve denotes the cross 
section for the nd breakup reaction. 



we can derive a relation between and 



Here, amplitudes (i = 0, 1, 2, . . . ) are defined as 



(AH) 



which satisfy 



Tin = 



(jp[o]| a \F® 



(Fl°]| A - T\- i+1 ~\ ' 



(A12) 



(A13) 



Calculations of a N-th order approximation start by 
regarding 



and thereby 



TW = (pio] 



(A14) 



(A15) 



Then, using Eqs. (|A11[) and (|A13|) backward, we calcu- 
late l*^ -1 !), |$[ w ^ 2 ]), successively until as 
the N-th order approximation for |$). 
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FIG. 14: (Color online) Differential cross sections of pd and nd breakup reactions for (a) the COL configuration, (b) the FSI 
configuration, (c) the SST configuration, and (d) the QFS configuration at Em = 13.0 MeV. The bold curves are for pd and 
thin curves for nd scattering. The dashed curves denote the calculation with the AV18 potential, and the solid curves those 
with the AV18+BR.660 potential. Experimental data are from Ref. [4l| (solid squares) for pd and Ref. [Z^] (open circles) and 
Ref. [53] (solid circles) for nd scattering. The arrows indicate the kinematical points that match the typical configurations. 



APPENDIX B: COULOMB FUNCTIONS 

In this appendix, we summarize formulae of functions 
related to spectator functions modified by the Coulomb 
potential. See Ref. [2l[ for details. 

Let Fi(r)(p),py) and G(,(r)(p),py) be the regular and 
irregular Coulomb functions that satisfy 
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ye(v(p),py) 



37T 
4m' 



(Bl) 

where ye(r)(p),py) is ether Fe(r](p),py) or Ge(r)(p),py), 
and 

2m e 2 



2 

A scattering state for the Coulomb potential — with 



energy j^p is written as 



F c (y;p, V (p)) = 4nJ2i% m *(Pm m (y) 



£,m 



xc t*,.MvV F t(v(p)>PV) (B3) 

py 



where er^ry) is the Coulomb phase shift, 
cr £ (ry) = ar g r(£+ 1 + vq). 



(B4) 



Analytical for of the Green's function, Eq. 
given by 



. 4m e^Wri^+^fc),^) 

GcAy,y,E P ) = -— 2 p — 



Fe(y(p),py<) 
py< 



(B5) 



where 11^(77, r) is defined as 

uf'fo.r) - e^'M (G e (i 1 ,r)±iF l ( V ,r)), (B6) 
giving the asymptotic form as 

u { ( ±] {ri,r) -> exp(±«(r-r/ln2r-^7r/2)). (B7) 

r — ►oo 

APPENDIX C: GREEN'S OPERATOR 

In this appendix, we first review two-body Green's op- 
erators, and then describe how to calculate Eq. (|3"3"|) for 
the case of k — 3 and E q > 0. 

We define Green's operators for the outgoing (+) and 
the incoming (— ) boundary conditions with a potential 
consisting of a short-range potential V s (x) and a long- 
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range Coulomb potential V c (x) as 



G 



(±) _ 



E q ±ie- T L [x) - V s (x) - V c (x) ' 



G 



(±) 

C,L 



E q ±iE-T L (x)-V c (x)' 



(CI) 



(C2) 



which satisfy resolvent relations 
r (±) _ r {±) r {±) v s r {±) _ r {±) , r {±) v s r {±) 

b i -^CL+^i v U C,L — «C,£ Lr C,L v U L ■ 

Two-body scattering wave functions corresponding to 
the outgoing and the incoming boundary conditions 
l^)^) satisfy the (partial- wave) Lippmann-Schwinger 
equations 

\4 ±) ) = \F L )+G { ±iv s \4 ±} ), (C4) 

where Fi,(^{q),qx) is a reduced Coulomb function de- 
fined by 

F Mq Uz) = FMq) ' qx) (C5) 



<y.f; 



with 



From the standing wave solution, the outgoing and the 
incoming solutions are obtained as 



1 



where /Cl is the scattering i^-matrix defined by 

777 A — 

= -q-^{F L \V\i, L ), 



(C13) 



(C14) 



which becomes tan 5 with a phase shift parameter 5. Us- 
ing the relations above, one obtains a relation between 
G^ ±) and VG L as 



G 



(±) 



1 



(C15) 



which reduces to Eq. (|C10|) if V(x) was 0, leading to 
^l{x) = FL{ r i{q),qx) and Kl = 0. 

Next, we discuss about asymptotic form of the Green's 
functions. 

Using the resolvent equation (|C3[) . the formal solu- 
tions, Eqs. (|C7|) and (|C12|) . and the asymptotic forms 
of the Coulomb Green's functions, which are obtained 
from their analytical forms, 



G 



(±) 



(C16) 



7(a) 



2h 2 q 



(C6) 



For later use, we define reduced Coulomb functions 
Gh{l{q),qx) and (7(g), qx) similar to Eq. (|C5jl . 

Using Eq. (|C3[) . we see that a formal solution of Eq. 
(IC4I) is written as 



ivi ±) ) = i^) + 4 ± v s iF L ). 



(C7) 



It is convenient to use the principal values of the two- 
body Green's operators defined as 



VG L =V 



1 



E q - T L (x) - V s (x) - V c (x) ' 



VG c ,l = V 



E q -T L {x)-V c {x)' 



(C8) 



(C9) 



As is Gq L , the analytical form of PG^ L is known and 
these operators are related as 



G { cl =, PGc,LTiq^\F L ){F L \ 



(CIO) 



A scattering wave function corresponding to VGc,l, 
namely standing wave solution \iPl) satisfies 

= \F L ) + VG C xV s \^l), (Cll) 
and a formal solution of this is given as 

\^l) = \F l )+VG l V s \F l ). (C12) 



VGc,L^q-^\G L )(F L \, 
n 



(C17) 



we obtain the asymptotic form of the Green's functions, 



G 



(±) 



-q^e^\u£W?\ 



VG L ^q-^\G L )(^ L \ 
n 



(C18) 



(C19) 



Finally, we describe how to calculate Eq. (f3"3"|) for k = 3, 
which we write simply as 

6{x) = (x\G ( + } \uj). (C20) 
Using Eq. (|C15[) . one can write 9(x) as 

— m — I — 
6( x ) = 6(x)-tq^ L (x)- — & L \u), (C21) 

where a new function 9{x) is defined by 

9(x) = (x\VG l \uj). (C22) 



From Eq. (|C19|) . the asymptotic form of 6(x) can be writ- 
ten as 

— TYl " — 

0(x) - q-^G L ( 7 (q),qx)(i; L \u>) (C23) 

In actual calculation, the function 9{x) is obtained by 
solving the ordinary differential equation 

[E q - T L {x) - V s {x) - V c {x)\ 0(x) = lo{x) (C24) 
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with the boundary condition 

6{x) cx 6 L (7(q),qx). (C25) 

2?— 'CO 

These relations give the asymptotic form of 9(x) as 

^(x)^ e+ -4 +) (7(.),^ )T -L- (- ff ™) (M.). 

(C26) 
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